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PREFACE 

The  causes  of  polar  auroras  and  electrojet  current  systems  have 
been  the  subject  of  much  conjecture  by  workers  in  the  field  of  geophysics. 
In  the  present  report,  the  author  attempts  to  relate  these  phenomena  to 
the  effects  of  solar  streams  on  the  regions  of  geomagnetically  trapped 
radiation . 

This  Memorandum,  which  extends  the  ideas  presented  in  RM-2753-NASA, 
Geomagnetic  Field  Distortion  by  a  Solar  Stream  as  a  Mechanism  for  the 
Production  of  Polar  Aurora  and  Electro jets,  represents  part  of  a 
continuing  study  of  the  properties  of  charged  particles  and  fields  in 
space,  and  was  prepared  for  publication  in  the  Journal  of  Geophysical 
Research. 

This  research  is  sponsored  by  the  National  Aeronautics  and  Space 
Administration  under  Contract  No.  NASr-2l(05)  and  monitored  by  the  NASA 
Office  of  Space  Sciences. 
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ABSTOACT 

A  mechanism  for  charge  separation  in  the  geomagnetically  trapped 
radiation  is  described  vhich  may  account  for  some  observed  phenomena 
associated  with  the  polar  aurora  and  the  electrojet  current  systems. 
Surfaces  of  constant  number  density  may  be  separated  from  surfaces  of 
constant  integral  invariant  within  the  trapped  radiation  as  a  result  of 
distortion  of  the  geomagnetic  field  by  solar  streams.  Drift  separation 
of  protons  and  electrons  will  follow,  and  for  irregular  distributions 
of  plasma  number  density,  electric  fields  will  arise.  A  direct  conse¬ 
quence  of  such  polarization  of  the  geomagnet ically  trapped  radiation 
will  be  the  polar-electrojet  current  systems. 

ftie  polar  aurora  arise  where  energetic  particles  are  discharged 
from  regions  of  excess  charge  within  the  geomagnetically  trapped 
radiation.  A  model  for  the  discharge  of  such  auroral  particles  iB 
discussed. 

An  interesting  feature  of  the  proposed  mechanism  is  that  the 
extreme  thinness  of  auroral  sheetB  appears  to  follow  as  a  natural 


consequence  of  charge  separation.  Uiis  is  shown  analytically  for  a 
simple  two-dimensional  model  of  a  trapped  plasma. 
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I.  INTRODUCTION  AND  SUMMARY 

In  a  letter  to  the  editor  of  the  J.  Geophys.  Research  (Kern,  1961), 
the  author  suggested  the  possibility  that  solar-stream  distortion  of  the 
geomagnetic  field  might  lead  to  a  situation  in  which  charge  separation 
occurs  in  trapped  radiation,  thus  giving  rise  to  polar-electro jet  current 
systems.  The  physical  basis  for  the  mechanism  suggested  is  adiabatic 
particle  motion  caused  by  eastward-  or  westward-directed  geomagnetic  field 
gradients.  Such  geomagnetic  field  gradients  could  be  expected  as  a 
result  of  the  interaction  of  the  geomagnetic  field,  containing  trapped 
plasma,  and  an  enveloping  ionized  solar  stream  of  comparatively  low- 
energy  particles. 

Charge  separation  had  been  proposed  earlier  by  Vestlne  (i960), 
Chamberlain.  Kern,  and  Vestlne  (i960),  and  Kern  and  VeBtlne  (1961)  to 
account  for  some  of  the  observed  features  of  polar  auroras  and  electrojet 
current  systems.  Hie  present  paper  discusses  the  mechanism  for  charge 
separation  outlined  in  the  author's  letter  (Kern,  1961)  and  discussed  by 
Chamberlain  (1961)  in  connection  with  his  auroral  theory.  In  particular, 
the  effects  of  increasing  electric  fields  are  studied.  It  is  found  that 
thin  auroral  sheets  can  be  accounted  for  by  the  proposed  charge  separation. 
Discharge  of  auroral  particles  from  regions  of  excess  charge  in  the 
trapped  plasma  is  also  discussed. 

Distortion  of  the  geomagnetic  field  by  ionized  solar  streams  has 
been  considered  by  many  authors  in  relation  to  the  occurrence  of  polar - 
electrojet  current  systems  and  aurora.  Hie  magnetic-storm  theories  of 
Chapman  and  Ferraro  (1931>  1932,  1933)  and  Martyn  (1951)  require  such 
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distortion.  Modification  of  the  geomagnetic  field  by  diamagnetic  trapped- 
particle  moments  and  by  current  systems  has  also  been  considered  by 
Dessler  and  Parker  (1959)/  Akasofu  (i960)  and  others.  Recent  progress  in 
the  understanding  of  relationships  between  plasmas  and  magnetic  fields 
has  led  to  several  pictures  of  the  form  of  a  solar- stream-distorted 
geomagnetic  field  (Johnson,  i960;  Piddlngton,  1959/  I960;  Parker,  i960). 

15ie  mechanism  described  in  the  author's  letter  (Kern,  1961)  is 
similar  in  certain  respects  to  that  described  by  Martyn  in  his  theory  of 
magnetic  storms  and  aurora  (Martyn,  1951)/  based  on  Chapman  and  Ferraro's 
ring-current  theory  of  magnetic  storms  (Chapman  and  Ferraro,  1931/  1932, 
1933).  In  the  present  paper,  however,  the  polarization  of  geomagnetically 
trapped  radiation  is  ascribed  to  the  displacement  of  particle  density 
distributions  with  respect  to  surfaces  of  constant  integral  invariant. 

Two  alternative  descriptions  of  the  mechanism  for  charge  separation 
in  a  geomagnetically  trapped  plasma  are  proposed  here.  First,  plasma 
density  irregularities  somehow  become  inclined  with  respect  to  surfaces 
of  constant  integral  invariant  I;  hence,  the  adiabatic  drift  of  protons 
westward  and  electrons  eastward  along  surfaces  of  constant  I  leads  to 
charge  separation.  Or,  second,  a  component  of  the  geomagnetic  field 
gradient  arising  from  distortion  of  the  geomagnetic  field  becomes 
oriented  parallel  to  surfaces  of  constant  plasma  density. 

Chamberlain  (1961)  has  studied  the  effects  of  charge  separation  in 
geomagnetically  trapped  radiation  in  relation  to  auroral  displays.  His 
treatment  is  based  on  conservation  of  the  total  energy  of  a  particle 
and  relates  change  of  electrostatic  potential  to  equal  and  opposite 
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changes  in  the  transverse  kinetic  energy  of  trapped  particles.  In  his 
model,  acceleration  of  particles  parallel  to  the  magnetic  field  Bf is 
associated  with  recovery  of  the  original  kinetic  energy  in  motion  along 
B.  The  treatment  is  applied  to  ray-type  auroral  structures.  Hie  present 
development  will  apply  to  auroral  features  of  the  sheet  type;  e.g.,  homo¬ 
geneous  arcs,  draperies,  and  bands. 

In  Section  II,  charge  separation  resulting  from  adiabatic  particle 
motions  in  the  geomagnetically  trapped  radiation  is  examined.  A  two- 
dimensional  model  is  analyzed,  and  it  is  found  that  the  extremely  thin 
sheets  characteristic  of  auroral  structures  may  be  a  natural  consequence 
of  such  charge  separation. 

In  Section  III,  the  effects  of  an  electric-field  component  parallel 
to  the  magnetic  fieldJ5  are  examined.  The  first  particles  discharged 
from  a  flux  tube  containing  excess  charge  will  be  those  with  small  magnetic 
moments.  As  discharge  continues,  the  magnetic -moment  distribution  of 
particles  within  a  flux  tube  will  tend  to  have  a  lower  bound.  This 
lower  bound  of  the  magnetic-moment  distribution  in  a  region  of  excess 
charge  is  used  as  a  basis  for  defining  a  model  electric-field  component 
parallel  to  J3.  A  model  electric  field  is  defined  which  has  Just  that 
magnitude  and  sense  which  permits  those  particles  of  one  sign  with 
lowest  magnetic  moment  to  move  freely  along  a  magnetic  field  line. 

The  motions  of  charged  atmospheric  particles  up  a  flux  tube  in 
the  presence  of  such  an  electric  field  are  examined.  Time  scales  for 
proton  neutralization  of  an  electron-charge  excess  are  found  to  be  long 
compared  to  electron  discharge  times.  Thus  electron  precipitation  can 
be  anticipated.  The  lifetime  of  a  proton-charge  excess  has  been  given 
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by  Chamberlain  (1961)  as  about  1  sec. 

Effects  of  the  model  electric  field  on  mirror  points  of  particles 
and  on  mirror-point  energies  are  computed.  Rates  of  change  of  the  model 
electric  field  are  computed  for  a  constant  flux  of  auroral  particles  into 
the  atmosphere  from  a  trapped  plasma.  The  time  scale  for  such  an  auroral 
flux  is  derived  and  is  al  :t  10  ^  sec  for  a  plasma  density  of  10  particles/cm^ 
and  a  flux  density  in  the  atmosphere  of  about  10^  particles/ cm2  sec.  The 
same  time  scale  applies  to  the  growth  of  the  electric  field  along  B.  One 
significant  result  of  the  analysis  is  that  particles  mirroring  at  a 
given  altitude  will  increase  in  energy  as  the  display  proceeds. 

The  relationships  between  charge  separation  in  a  geomagnet ically 
trapped  plasma  and  ionospheric  current  systems  are  discussed  in  Section 
IV.  An  estimate  is  made  of  the  distortion  of  the  geomagnetic  field 
necessary  to  produce  a  typical  westward-directed  electrojet.  It  is 
concluded  that  moderate  distortions  would  suffice  in  the  presence  of 
plasma-energy  densities  that  are  clearly  less  than  the  geomagnetic-field 
energy  density  in  regions  connected  with  the  auroral  zone. 

Charge  separation  in  the  geomagnet ically  trapped  radiation  has 
implications  for  auroral  morphology,  and  these  are  discussed  in  Section  V. 
Correlations  with  observations  are  discussed,  and  a  possible  means  of 
generating  horseshoe  auroral  arcs  is  outlined.  Section  VI  deals  with  a 
possible  explanation  of  the  tendencies  of  polar  auroras  and  electrojet 
current  systems  to  recur. 
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II.  CHARGE  SEPARATION  AND  POLARIZATION 
IN  GEOMAGNETICAL!*  TRAPPED  RADIATION 

In  the  geomagneti cally  trapped  radiation,  there  are  two  distinct 
types  of  adiabatic  drift  characterized  by  whether  protons  and  electrons 
drift  in  the  same  or  in  opposite  directions.  JgxB, drift  or  any  drift 
caused  by  forces  that  depend  on  the  sign  of  electric  charge  is  in  the 
same  direction  for  both  protons  and  electrons.  Inertial  forces  or  those 
dependent  on  magnetic  moment  will  lead  to  adiabatic  drift  of  protons  and 
electrons  in  opposite  directions. 

The  familiar  westward  drift  of  trapped  protons  and  eastward  drift 
of  trapped  electrons  due  to  the  curvature  and  gradient  of  the  geomagnetic 
field  is  given  by 

V.B  XB  2 

v  - - - -  W(l  +  cos  a)  (1) 

e  B3 

where  V.  B  is  the  component  ofVB  transverse  to  B,  e  is  the  charge  of  the 

X  ^ 

particle  in  emu,  W  is  the  kinetic  energy  of  the  particle  in  ergs,  and  a 
is  the  pitch  angle  of  particle  motion  with  respect  tojg,  (Spitzer,  1956). 

In  a  plasma  with  an  irregular  number  density,  this  kind  of  drift  can 
lead  to  charge  separation  with  associated  electric  fields.  Figure  1 
shows  the  result  of  this  kind  of  charge  separation  in  a  bounded  sheet  of 
plasma. 

In  applying  the  idea  of  charge  separation  to  the  geomagneti cally 
trapped  radiation,  distortion  of  the  field-plasma  geometry  must  be 
introduced  in  such  a  way  that  the  drift  velocity  is  not  simply  parallel 
to  surfaces  of  constant  number  density.  Figure  2  shows  a  configuration 
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in  which  a  component  of  (denoted  VjpB)  is  parallel  to  surfaces  of 

constant  number  density.  This  component  of  leads  to  charge  separation 

and  an  electric  field  directed  from  the  separated  proton  sheet  toward  the 
electron  sheet. 

It  can  be  seen  that  if  the  transverse  electric  field  arising  from 
charge  separation  is  increasing,  the  JSxB^ drift  of  both  protons  and 
electrons  will  be  accelerated.  Chandrasekhar  (i960)  shows  that  the 
adiabatic  character  of  the  particle  motion  is  preserved  for  this  situation 
if  such  a  polarization  drift  velocity  v  is  introduced  that 

O  c>E 

=  (m/eB  >  3?  <2) 

where  m  is  the  particle  mass,  e  is  the  electronic  charge  in  emu,  B  is 

•dE 

the  magnetic  field,  and  is  the  rate  of  change  of  the  electric  field 

in  emu/sec.  Uiis  drift  exactly  compensates  for  the  increase  in  kinetic 

energy  arising  from  acceleration  of  the  JSxB^  drift  by  a  simultaneous 

decrease  in  electrostatic  potential  energy.  Note  that  this  drift  is 

independent  of  particle  energy  and  is  much  faster  for  protons  than 

electrons,  because  of  the  larger  proton  mass.  For  a  trapped  plasma  of 

the  kind  discussed  here,  the  polarization  drift  v^  opposes  the  magnetic 

drift  v  .  A  two-dimensional  plasma  model  will  therefore  be  described  for 
B 

which  the  electric- field  variation  and  its  effects  can  be  calculated. 

Consider  a  two-dimensional  plasma,  such  as  that  indicated  in  Fig.  1, 
for  which  surfaces  of  constant  plasma  density  contain  magnetic  field 
lines.  Suppose  a  magnetic  field  gradient  ^B  exists  transverse  to  the 
magnetic  field  B.  Hie  drift  motions  of  protons  and  electrons  in  the 
plasma  that  are  due  to  this  gradient  will  be  supposed  to  have  a  component 
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perpendicular  to  surfaces  of  constant  plasma  density.  In  order  to 

simplify  the  calculations,  it  will  be  assumed  that  the  proton  drift  is 

much  greater  than  that  of  the  electrons.  This  would  follow  from  (l)  if 

the  electrons’  kinetic  energies  were  very  much  less  than  those  of  the 

protons.  If  this  is  the  case,  motion  of  a  proton  can  be  specified  locally 

by  the  drift  velocity  v,  =  v_  -  v  ,  where  vl  and  v  are  given  by  Eqs.  (l) 

and  (2),  respectively.  For  the  two-dimensional  geometry  considered 

here,  Gauss*  theorem  can  be  applied  to  obtain  the  electric  field  at  any 

point  inside  the  plasma.  Here  the  velocity  v_  will  be  considered  to 

be  the  component  of  the  gradient  drift  velocity  that  is  perpendicular  to 

surfaces  of  constant  plasma  density.  Such  drift  would  result  from  a 

component  of  parallel  to  the  surfaces  of  constant  plasma  density. 

Right-handed  coordinate  axes  x,  y,  and  z  are  taken  so  that  the  z 

axis  is  in  the  direction  of  the  magnetic  field  B.  Surfaces  of  constant 

plasma  density  are  taken  parallel  to  the  x-z  plane,  and  the  y-axis  is  in 

the  direction  of  the  plasma  density  gradient  n(y).  !Hie  plasma  number 

density  is  taken  as  varying  from  zero  to  some  value  n(y)  as  y  goes  from 

-®  to  y.  A  gradient  inj^  in  the  -x  direction  will  lead  to  proton  drift 

v  in  the  -y  direction.  The  polarization  drift  v  of  a  proton  will  be 
B  P 

in  the  +y  direction;  that  is,  v^  is  in  the  direction  of  increasing  n(y). 

The  total  drift  velocity  v,  =  v_  -  v  will  be  in  the  -y  direction  if 

d  B  p 

v  >  v  ,  and  in  the  +y  direction  if  v  >  v  .  Ihe  drift  v  is  in  such  a 
Bp  p  B  D 

direction  that  protons  will  tend  to  separate  out  of  the  neutral  plasma, 
as  shown  in  Fig.  1.  On  the  opposite  side  of  the  neutral-plasma  feature, 
electrons  will  separate  out.  An  electric  field  E  will  appear  in  the  +y 


direction. 
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Now  the  rate  of  change  of  space  charge  per  unit  volume  resulting 
from  the  drift  velocity  is  given  by 


!l  -  .1”  t4<w)  «  (3) 

VO 

where  v^(w)  is  the  total  drift  velocity,  which  is  a  function  of  proton 
kinetic  energy  W,  and  F(e,y)  is  a  distribution  function  for  total  proton 
energy  e,  so  that 

p  00 

I  F(e,y)  de  =  n(y). 

Bie  total  energy  e  is  invariant  for  a  given  proton.  Hie  kinetic  energy 
W  will  vary  with  the  development  of  electrostatic  potential  V,  so  that 
W  =  €  -  eV. 

2 

Now  for  the  geometry  selected  ^  •  E  =  4rrc  p  (in  emu)  reduces  to 


and  hence 


dE 


hire 


2 

P 


From  Eq.  (3)  it  follows  that 

VB(W)  '  Vp]  d€  (k} 

Substituting  from  Eqs.  (l)  and  (2)  and  noting  that  W  =  €  -  eV,  Eq.  (4) 


can  be  written 
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$=Sf  -  ^c2P|  (e  -  ev)  d.  -  iffifs!  | 


.  .-4TJ1 

where,  from  Eq.  (l),  p  =  v^/W  =  constant.  Note  that  ^  and  V  Eire  not 

explicitly  functions  of  e.  Also 

roo 


while,  if  e  does  not  depend  on  y. 


d€  =  ?  |s 

dy  dy 


where  I  is  the  average  total  particle  energy.  Hence,  Eq.  (4)  can  be 
written 


2  2 
d  E  , _ 2n  „s  dn  4innc  dn  dE 

■  W  6  <e  -  'v>  ■  -2-  3?  3t 

D 


d  E  .  4Trmc  dn  dE  ,  2„  „x  dn 

3y5t  *  -T-  3t  ‘  13  (E  -  eV) 


In  order  to  solve  this  differential  equation,  the  condition  that 
eV  «  e  will  be  imposed.  Ihis  will  be  valid  initially.  For  later 
times,  the  growth  of  the  electric  field  will  be  retarded  as  eV  becomes 
significant  with  respect  to  the  mean  particle  energy.  If  eV  «  e  the 
right  side  of  (6)  can  be  written  47rec  v^  ,  where  v^  is  the  mean 
drift  velocity  due  to  the  magnetic  field  gradient.  Since  the  energy 


distribution  is  taken  as  independent  of  y,  v  is  independent  of  y. 

B 
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Equation  (6),  then,  has  the  solution 


¥ 


I  2  -  Bn 
kirec  vb37 


e^ 


2  2 

where  y  =  l+Trnc  /B  .  This  equation  can  now  be  integrated  with  respect  to 

y  from  some  point  outside  the  region  of  the  plasma  density  feature 

<3e 

(where  both  n  and  are  zero)  to  y.  The  result  is 


Be 


e7n  =  (4frec 


2 


vB/r)(eyn 


1) 


Hence  the  rate  of  change  of  the  electric  field  is  given  by 


Be 

3t  = 


Itirec^v, 


B 


) 


(7) 


where  e  is  the  proton  charge  in  emu,  c  is  the  velocity  of  light,  v_  is 
the  mean  drift  velocity  for  protons  in  the  negative  y-direction 
resulting  from  a  transverse  magnetic- field  gradient,  n  is  the  neutral 
plasma  number  density,  and  y  =  kirmc  /b  ,  where  m  is  the  proton  mass  and 
B  is  the  value  of  the  magnetic  field.  This  expression  is  seen  as  valid 
initially.  As  mentioned  above,  the  development  of  an  electrostatic 
potential  field,  V,  will  modify  the  drift  motion,  effectively  reducing 

Be 

the  mean  drift  velocity  v^.  This  will  lead  to  a  decrease  in  .  The 
electric  field  does  not,  therefore,  increase  without  limit. 

The  rate  of  change  of  charge  density  p  associated  with  this 

Be  2 

distribution  of  electric  field  can  be  computed  from  ^  *=  krrc  p,  which 


B2E 

Bydt 


gives 
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or,  from  (7) 


5t  =  evB  6 


•yn  dn 
<3y 


(8) 


This  expression  gives  the  time  variation  of  space  charge  p  subject  to 
the  same  restrictions  as  discussed  for  as  given  by  (7).  The  space 
charge  p  is  initially  distributed  very  near  the  edges  of  any  assumed 
plasma  density  feature.  The  extreme  thinness  of  such  a  space  charge 
distribution  can  be  emphasized  by  considering  an  example. 

For  a  plasma  of  number  density  n,  the  value  of  y  for  protons  is 

-2  2 

about  I.87  x  10  B  ,  where  B  is  the  intensity  of  the  magnetic  field. 
Thus,  for  values  of  B  applicable  to  the  geomagnetic  field  connected  with 
the  auroral  zone  of,  say  4  x  10  J  gauss,  y  ~  10  .  It  is  apparent  from 
(7)  that  would  be  nearly  zero  on  the  interior  of  a  plasma  density 
feature  where  n  becomes  significantly  different  from  zero.  Talcing  a 
single  plasma  density  distribution  of  the  form  n(y)  =  a  sin  (=r)y, 
where  a  is  the  maximum  n  and  L  is  the  width  of  the  distribution,  leads 

to  the  result  that  the  maximum  of  ^  occurs  within  a  distance  of 

2  2 

about  2L i/r  on  the  edge  of  the  distribution,  where  y  =  l+Trmc  a/B  . 

'  'max  ^  'max  ' 

Fbr  b  /v  x  10  ^  gauss  and  n^lO  particles/cm^,  the  maximum  of 
would  occur  within  a  distance  of  about  L/lO^  from  the  edge  of  the 
distribution.  The  regions  of  excess  charge  are  therefore  very  thin  and 
appear  very  near  minima  in  the  neutral  plasma  density  for  the  model 
analyzed  here.  The  same  conclusion  can  be  drawn  for  the  electron 
excess  corresponding  to  the  proton  excess  dealt  with  here.  The  constant 
7  derived  from  the  proton  motion  must  also  be  applied  to  the  formation 
of  an  electron  charge  excess,  since  polarization  drift  of  protons  will 
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affect  the  formation  of  an  electron  sheet. 

In  order  to  fully  illustrate  this  result,  the  thin-sheet  character 
of  the  excess-charge  sheath  can  he  derived  from  consideration  of  the 
electric  field.  Die  excess  charge  is  the  source  of  the  electric  field. 
Hence,  0.99  of  the  excess-charge  must  appear  between  the  edge  of  the 
distribution  and  a  point  inside  the  distribution  where  the  time  rate  of 
change  of  E  is  0.99  the  maximum  For  the  sample  plasma  number 

p  77” 

density  distribution  n(y)  =  a  sin  =-  y  considered  here,  a  sheath  thickness 

b  can  be  defined  as  that  distance  inside  the  distribution  which  includes 

0.99  of  the  excess  charge.  Dien  fran  (7),  (^)  occurs  for  n(^)  =  a. 

For  =  0.99  (||)  ,  1  -  e-711^  =  0.99  (l  -  e"7a).  In  general,  for 

max 

a  >  10-  particles/cm3,  e~7&  «  1.  Hence  e~7n  =  1  -  0.99  =  0.01,  and 
n(b)  =  (1/7)  in  100.  Since  7~10^,  n(b)  =  4.6  x  10-^  particle s/cm^. 

For  small  n(b)/a,  sin  Trb/L^TTb/L  =  >/n(b)/a,  hence  (b/L)  =  (1/7 r)  \/n(b)/a. 

For  a  =  10  particles/cm^,  this  gives  b/L^  f  x  10  Table  1  gives  a 
resume  of  b/L (as  defined  here  for  the  plasma  density  distribution 
n(y)  =  a  sin  (tt/l)  y)for  different  values  of  the  maximum  number  density 
a.  Also  shown  is  the  thickness  b  in  km  for  L  =  1000  km.  Die  values  of 
b/L  shown  here  would  be  considered  applicable  to  the  equatorial  plane. 
Projected  along  B  into  the  ionosphere  in  the  auroral  zone,  sheets  of 
thickness  b  would  be  reduced  in  thickness  by  a  factor  of  about  30*  Dius 
a  trapped  plasma  feature  1000  km  in  thickness  with  a  peak  density  of 
1  particle/cm^  in  the  equatorial  plane  might  lead  to  an  auroral  arc  with 
a  thickness  of  2/3  km.  Die  incident  pnrticle  flux  in  such  an  aurora 
might  typically  be  10  particles/cm  sec.  Diis  would  require  a  mean 
rate  of  charge  separation  of  10  charged  particles/sec  in  a  magnetic 
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Table  1 
RESUME  OF  b/L 

Ratio  b/L  of  Excess  Charge  Sheath  Thickness  b  to  Plasma  Density  Perturbation 
Thickness  L  for  Plasma  Density  n(y)  =  a  sin  (7r/L)y,  Where  a  is  Peak 

Plasma  Density 


Peak  Density  a 
Particles/ cm^ 

Ratio  of  Sheath  Thickness  to 
Perturbation  Thickness  b/L 

Sheath  Thickness  b 
for  L  =  1000  km 
Plasma  Density 
Perturbation,  km 

0.01 

0.22 

220 

0.1 

0.068 

68 

1 

0.022 

22 

10 

0.007 

7 

100 

0.002 

2 
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flux  tube  (connected  with  the  auroral  zone)  that  has  a  volume  of  about 
10  cm  .  From  Eq.  (8),  such  a  flux  could  be  supplied  by  protons 

drifting  with  a  mean  velocity  v  of  about  100  m/sec  and  distributed  so 

cLn  "3 

that  gradients  ^  occur  of  the  order  a/L  or  10  particles/cnr3/ 1000  km, 

In 

provided  ^  is  finite  very  close  to  the  zeros  in  n(y). 

A  background  plasma  would  suppress  ^  ,  as  can  be  seen  from  (8). 
Consider  n(y)  =  n  +  n  (y),  where  n  is  constant  and  n  (y)  is  the 
perturbation  plasma  density.  If  nQ  >  0,  (8)  can  be  written 

=  ev  e  0  e  P  This  rate  of  change  of  p  is  obviously  less 

dy  ~7n0 

than  that  which  would  hold  for  n^  alone  by  the  factor  e  resulting 

from  the  background  plasma.  This  effect  will  serve  to  reduce  the 
charge  separation  and  hence  any  resulting  auroral  particle  fluxes.  It 
can  therefore  be  concluded  that  for  auroral  displays  to  be  a  result  of 
charge  separation  in  the  magnetosphere,  the  background  plasma  must  be 
discontinuous.  If  structure  in  the  plasma  density  In  the  magnetosphere 
is  due  to  plasma  instabilities,  the  background  plasma  is  likely  to  have 
a  structure  at  least  as  complex  as  that  inferred  for  higher  energy 
particles  (Rosen,  Farley,  and  Sonett,  i960). 

One  inference  can  be  drawn  from  this  admittedly  oversimplified 
treatment  of  charge  separation  In  a  two-dimensional  plasma.  Ihat  is, 
if  the  electric  fields  resulting  from  charge  separation  are  increasing, 
there  will  be  a  tendency  for  excess  charge  to  appear  in  fairly  restricted 
regions  near  the  outer  boundaries  of  plasma  density  perturbations.  Thus 
it  seems  possible  to  account  for  the  extreme  thinness  of  auroral  sheets 
within  the  framework  of  an  auroral  theory  based  on  charge  separation 
resulting  from  adiabatic  particle  motion. 
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It  will  be  seen  in  Section  III  that  increasing  electric  fields  are 
also  required  to  account  for  continuing  discharge  of  auroral  particles 
from  a  geomagnet ically  trapped  plasma. 


% 
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III.  glE  LOWERING  OF  THE  MIRROR  POINTS  OF  TRAPPED  PARTICLES 

When  charges  separate,  as  described  above,  reflection  of  the 
trapped  particles  occurs  at  lower  altitudes,  provided  that  electric 
fields  arise  parallel  to  the  magnetic  field  lines.  As  will  be  suggested 
in  Section  IV,  electrojett.  my  result  when  thermal  particles  neutralize 
excess  charge  of  opposite  sign,  but  fluxes  of  more  energetic  particles 
Into  the  atmosphere  are  required  to  form  auroral  displays.  As  pointed 
out  by  Kern  and  Vestine  (lS>6l)  and  Chamberlain  (1961),  many  of  the 
features  of  auroral  morphology  can  be  accounted  for,  if  excess  charge 
distributions  can  be  maintained  in  the  trapped  radiation. 

Chamberlain  (1961)  discusses  the  lowering  of  mirror  points  in 
regions  of  excess  charge.  The  drift  of  particles  antiparallel  to  the 
direction  of  the  transverse  electric  force  on  particles  leads  to  the 

conversion  of  transverse  kinetic  energy  into  electrostatic  potential 

•  • 

energy  at  the  rate  W,  =  -  eE  '  v.  =  -  eV,  where  eE  Is  the  electric 

J.  ****&  AV 

force  on  the  particle,  v  is  its  drift  velocity,  and  where  eV  is  the 
rate  of  change  of  electrostatic  potential  energy  of  the  particle.  Thus 
the  pitch  angle  of  a  particle,  relative  to  the  magnetic  field,  changes 
as  its  electrostatic  potential  changes.  Discharge  of  a  particle  from  a 
region  of  high  electrostatic  potential  (a  region  of  excess  charge)  along 
magnetic  field  lines  leads  to  a  recovery  of  this  kinetic  energy.  Hie 
system  considered  by  Chamberlain  is  a  flux  tube,  containing  excess 
charge,  of  about  the  same  dimensions  as  the  Larmor  radius  of  the  trapped 
particles.  In  this  configuration,  the  transverse  electric  field  of  the 
excess  charge  varies  in  direction  during  a  single  orbital  cycle.  Hie 
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magnet  ic  moments  of  particles  are  not  invariant  for  an  electric  field 

WA 

of  this  geometry,  although  =—  still  is. 

For  the  two-dimensional  geometry  considered  here  (Fig.  2),  the 
transverse  electric  field  is  essentially  constant  in  direction  and 
magnitude  during  a  single  orbital  cycle.  In  this  case,  although  the 
transverse  kinetic  energy  changes  while  the  particles  drift,  the 
particles'  magnetic  moment  is  invariant  ( Hellwig ,  1955)*  In  fact,  it 

can  be  shown  that  the  time  rate  of  change  of  the  transverse 

•  • 

kinetic  energy  is  given  by  =  |iB,  where  p  is  the  particle's  magnetic 

moment,  and  B  the  time  rate  of  change  of  B  due  to  Ep^  drift  motion  in 
a  nonuniform  magnetic  field  B.  Diis  is  equivalent  to  a  betatron 

MV 

deceleration  of  the  transverse  motion. 

In  regions  of  excess  charge,  accelerations  of  particle  motion 
parallel  to  may  occur.  In  order  to  relate  auroral  displays  to  particle 
fluxes  along  J3,  a  model  electric  field  will  be  defined.  Biis  electric 
field  is  assumed  to  arise  frcm  excess  charge  resulting  from  the  charge 
separation  mechanism  as  described  in  Section  II.  The  effects  of  such  an 
electric  field  on  a  trapped  plasma  will  be  examined. 

The  effect  of  a  model  electric  field  component  E2  parallel  to  the 
magnetic  field  B  can  be  examined  using  the  equations  of  motion  for  a 
single  particle.  In  general,  there  will  also  be  a  transverse  component 
E^j  this  will  be  neglected  here.  The  acceleration  of  a  particle  V 
parallel  to  B  will  be  given  by 

*  -  (l/m)(-n  +  eE2)  (9) 

where  m  is  the  masB  of  the  particle,  p  the  magnetic  moment,  the 
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component  of  the  gradient  of  B  parallel  to  B,  e  the  charge  of  the  particle, 
and  E  the  component  of  the  electric  field  parallel  to  B.  Here  z  is  taken 

d.  ^ 

as  a  displacement  from  the  equatorial  plane  along  and  the  electric 
force  eEg  is  taken  in  the  same  direction.  Particles  of  opposite  charge 
can  he  treated  by  taking  the  opposite  sign  for  e,  along  with  the  appro¬ 
priate  mass  and  magnetic  moment.  Initiation  of  charge  separation  will 
lead  to  the  discharge  along  B  of  the  particles  with  the  smallest  magnetic 

/W 

moments  from  the  region  of  excess  charge.  As  such  particles  discharge, 
those  remaining  will  have  larger  moments.  At  any  given  time,  it  seems 
plausible  to  define  a  lower  bound  for  the  moments  of  the  remaining 
particles.  The  magnitude  of  this  smallest  magnetic  moment  ^  will 
evidently  increase  as  discharge  continues. 

The  field  required  to  permit  free  motion  of  particles  along is 
given  by  equating  z  to  0  in  (8)  for  particles  with  this  smallest 
moment  \i Wiis  is  Just 

E2=(pye)V2B  (10) 

This  electric  field  will  produce  a  flux  <J  of  particles  with  magnetic 
moments  \x ^  into  the  atmosphere.  The  flux  $  must  correspond  to  the 
total  drift  flux  of  excess  charge  into  the  sides  of  the  region. 

The  latitude  oscillation  time  for  a  particle  with  magnetic  moment 

U  „  can  be  computed  for  a  given  magnetic  field  line  and  mirror  field  B 

Hi  m 

in  the  absence  of  electric  fields.  Ihe  acceleration  of  such  a  particle 
travelling  from  the  mirror  point  to  the  equatorial  plane  is  given  by 
where  V2 B  is  "the  parallel  component  of  the  magnetic  field 
gradient  and  m  is  the  particle  mass.  Consider  now  a  thermal  particle 
of  opposite  sign  travelling  fYom  the  same  mirror  point  to  the  equatorial 
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plane  in  the  model  electric  field  E^.  The  acceleration  of  this  particle 
will  be  given  by  eE^/M,  where  e  is  the  charge  and  M  the  mass  of  the 
thermal  particle,  and  the  model  electric  field  is  given  by  (10).  It 


is  evident  that  a  thermal  particle  of  opposite  charge  undergoes  an 
acceleration  m/M  times  that  undergone  by  a  particle  with  moment 
mirroring  in  the  atmosphere  in  the  absence  of  electric  fields.  Hence 
such  a  thermal  particle  has  a  travel  time  from  the  atmosphere  to  the 
equatorial  plane  which  is  \j M/m  times  the  travel  time  of  a  particle  of 
opposite  sign  oscillating  between  mirror  points  in  the  absence  of  an 
electric  field. 

A  model  electric  field  sufficient  to  cause  6  Kev  electrons  with 

—8 

moment  n  =  2  x  10  emu  to  discharge  along  B  into  the  atmosphere  will 
therefore  accelerate  atmospheric  protons  to  the  equatorial  plane  in 

/l840  x,  where  t  is  the  travel  time  of  an  electron  with  moment 

-8 

H  =  2  x  10  emu  which  mirrors  in  the  atmosphere  in  the  absence  of  an 
electric  field.  Biis  time  t  is  of  the  order  of  1  sec.  Time  scales  for 
proton  neutralization  of  electron  charge  excesses  are  therefore  large, 
allowing  the  development  of  electron  excesses  that  can  discharge  to 
generate  auroral  displays.  Electrons  from  the  atmosphere  can  quickly 
neutralize  proton-charge  excesses.  Chamberlain  (1961)  estimates  lifetimes 
of  the  order  of  1  sec  for  proton  charge  excesses. 

Tfre  above  argument  indicates  that  transient  discharges  of  electrons 
can  be  expected.  Auroral  displays,  however,  can  persist  for  many 
minutes,  and  in  some  cases  hours.  It  is  therefore  important  to  estimate 
possible  ion  fluxes  upward  along  B  into  an  auroral  flux  tube.  !Rie  model 
electric  field  E^  will  lead  to  such  upward  motion  of  ions.  In  the 
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conducting  ionosphere,  the  ion  current  density  along  can  be  represented 
by  ^  =  a0i^2  =  ^ie  /m±v±)^2>  wbere  o  ^  =  N^e  /m  is  the  ion  conduc¬ 
tivity  along  B,  N.  is  the  local  number  density,  e  is  the  charge,  m.  is 
the  mass,  and  is  the  collision  frequency  of  the  ions.  If  the  ions 
considered  are  protens,  the  charge  is  the  electron  charge.  Now  the  net 
ion  flux  upward  is  limited  by  the  conductivity  along  B  near  the  bottom 
of  an  auroral  flux  tube.  Ion  motions  in  the  outer  parts  of  such  a  flux 
tube  will  act  to  redistribute  excess  charge  resulting  from  charge 
separation.  Charge  neutralization,  however,  will  occur  only  by  the 
introduction  of  a  net  ion  flux  into  the  tube.  Except  for  the  regions  in 
the  ionosphere  where  the  conductivity  transverse  to  B  is  significant, 
a  net  ion  flux  can  only  be  introduced  into  an  auroral  flux  tube  by 
conduction  along  B.  The  value  of  the  net  ion  flux  into  an  auroral  flux 
tube  can  be  obtained  if  the  electric  field  E2  along  is  specified  at, 
say,  the  Y^  maximum  in  the  ionosphere.  If  this  ion  flux  is  less  than  the 
rate  at  which  excess  charge  is  appearing  in  the  flux  tube,  the  electric 
field  can  increase  and  energetic  electrons  can  be  discharged.  The  ion 
conductivity  a  .  along  B  in  the  F0  region  of  the  ionosphere  at  an 

Ol  am  d 
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altitude  of  300  tan  is  of  the  order  of  10  emu.  The  electric  field  Eg 
that  will  discharge  6-Kev  electrons  to  this  altitude  is  about  2  x  10J  emu. 

Jbr  this  electric  field,  the  ion  current  density  is  therefore 

•12  3  2  2 

10"  •  2  x  10J  abamp/cm  =  2  x  10"'*  abamp/cm  .  This  corresponds  to  an 

11  2 

ion  flux  density  of  10  protons/cm  sec  upward  along  B.  This  net 
upward  ion  flux  is  very  nearly  the  same  as  the  flux  density  of  6-Kev 
electrons  observed  by  Mcllwain  (i960)  in  a  bright  auroral  arc.  Hence, 


charge  separation  rates  inferred  from  auroral  electron  fluxes  can  be 
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great  enough  to  equal  or  exceed  ion  charge  transport  upward  along  B.  If 
this  is  the  case,  the  electric  field  E^  along  B  can  increase  with  time 
and  continuous  discharge  of  trapped  electrons  is  possible. 

A  fictitious  magnetic  moment  p'  can  be  defined  for  a  particle.  Let 
p'  =  p  t  p^,  where  p  is  the  actual  particle  moment  (which  is  invariant), 
p^  is  as  defined  earlier.  The  plus  sign  is  taken  for  particles  with  the 
opposite  charge  of  the  excess  in  the  region,  the  minus  sign  is  taken  for 
particles  with  the  same  charge  as  that  of  the  excess.  Then,  the 

acceleration  parallel  to 3^  is  given  by  z  =  -  (l/m)  p'  VgB.  It  follows 

•  • 

that  z  *=  -  (l/m)(p  ±  p^)  ^ZjB.  This  is  the  equation  of  motion  for 
particles  moving  in  the  combined  electric  and  magnetic  fields.  Hence, 
particle  motion  is  modified  by  the  electric  field  Eg  in  the  Bame  way  as 
if  the  particle's  magnetic  moment  were  changed  and  no  electric  field 
were  present.  For  electrons  moving  in  a  flux  tube  containing  an  electron 
excess,  the  fictitious  magnetic  moment  of  an  electron  is  p'  ■  p  -  p^, 
where  p^  is  the  lowest  electron  moment  present,  and  p  is  the  real 
magnetic  moment.  Nete  that  for  a  constant  electric  field,  p^  is  constant, 
and  hence  p'  iB  invariant. 

(Jhanges  in  the  mirror  fields  of  particles  can  be  examined  using 
the  fictitious  magnetic  moment.  Die  kinetic  energy  of  a  particle  as  it 
crosses  the  equatorial  plane  is  modified  by  transverse  drift  in  the 
nonuniform  electric  and  magnetic  fields.  If  the  magnetic  moment  p  is 
constant,  the  kinetic  energy  at  a  new  equatorial  magnetic  field  can 
be  given  by  -  Wq  -  p(BQ  -  B^),  where  Wq  is  the  particle's  kinetic 
energy  prior  to  charge  separation,  and  Bq  is  the  initial  value  of  the 
magnetic  field  where  the  particle  crosses  the  equatorial  plane.  Applying 
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the  idea  of  the  fictitious  magnetic  moment,  note  that  the  new  mirror 

field  B’  in  the  presence  of  the  electric  field  E~  is  given  by  W’/n** 
m  d  o 

Also  note  that  W  =  piB  ,  where  B  is  the  original  mirror  field  of  the 
o  nr  m 

particle.  Thus  B^  =  (p/^')(B  -  Bq  +  B^),  where  all  of  the  quantities 
have  been  previously  defined.  It  is  evident  that  if  and  Bq  are 
specified,  the  new  mirror  field  B^  can  be  calculated.  For  the  case 
in  which  Bq  -  B^  approaches  zero,  or  is,  at  least,  small  compared  to  B^, 

B* 'S’  (u/fi*  )  B  (11) 

m  '  m 

Since  pi*  =  pi  the  new  mirror  field  is  given  by  B*  =  (l  1  p,/p)  1  B  , 

ju  m  Yj  ni 

where  the  minus  applies  to  particles  with  the  same  charge  as  the  excess 

charge  giving  rise  to  the  electric  field  along  ^B,  and  the  plus  applies 

to  particles  of  opposite  charge. 

The  mirror  point  energies  of  particles  can  be  computed  using  the 
invariance  of  the  real  magnetic  moment  pi.  The  kinetic  energy  of  a 
particle  at  the  new  mirror  field  is  given  by  W  =  piB^.  This  is  just 
v  Wq,  where  Wq  is  the  original  kinetic  energy,  and  B^  is  the 
original  mirror  field. 

The  auroral  flux  produced  by  an  electric  field  Eg,  as  defined 
above,  can  be  computed.  The  time  constants  for  particle  discharge  to 
the  atmosphere  can  thus  be  evaluated.  Consider  a  given  initial  dis¬ 
tribution  of  trapped-particle  energies  within  a  magnetic  flux  tube 
linking  the  auroral  zone. 

In  the  following  treatment,  it  will  be  assumed  that  the  particles 
trapped  within  a  given  magnetic  flux  tube  are  not  resupplied  by 
adiabatic  drift.  The  source  of  electrostatic  energy  is  assumed  to  be  an 


% 
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excess  of  trapped  particles  of  one  charge  appearing  in  the  flux  tube  as 
a  result  of  charge  separation.  An  auroral  flux  will  thus  be  due  to 
discharge  of  particles  initially  contained  in  the  flux  tube.  This  can 
be  seen  to  occur  for  the  case  in  which  more  energetic  particles  of  one 
kind  (say,  protons)  drift  out  of  a  flux  tube,  leaving  an  excess  of 
particles  of  opposite  charge. 

The  polarization  component  of  adiabatic  drift  motion  which  results 
from  increasing  electric  fields  transverse  to  B  described  in  Section  II 
would  inhibit  the  drift  of  low  energy . particles  into  a  region  of 
separated  charge.  Thus  it  might  be  expected  that  any  drift  supply  of 
particles  to  regions  of  excess  charge  would  be  mainly  higher  energy 
particles.  If  this  is  the  case,  the  process  described  here  would  give 
the  auroral  flux  from  a  background  plasma  in  a  region  of  separating 
charge . 

Fbr  an  isotropic  distribution  of  particle  velocities  at  the 
equatorial  plane,  the  differential  particle-density  distribution  can 
be  given  by 

dn  *>  n  g(w)  dW  d  cos  a  (12) 

where  g(w)  is  the  distribution  function  for  particle  kinetic-energies, 

a  is  the  pitch  angle  and  W  the  kinetic-energy  for  a  given  particle,  and 

n  is  the  total  particle  density.  Now  W  =  nB^,  and  dW  -  B  dji,  for  a 

given  pitch  angle  a,  where  u  is  the  particle  magnetic  moment  and  Bm  is 

the  mirror  field.  Also  sin  a  =  B  /b  ,  whehe  B  is  the  equatorial  field. 

o'  m  o 

Hence,  d  cos  a  =  J  (B  /B2)(l  -  B  /B  )'^2  dB  .  It  follows  from  (12) 
d  or  m  o  m  m 


that 
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dn  =  (n/2)  g(^Bm)(Bo/lm)(l  -  B^Bj'1/2  dBm  d*i  (13) 

The  total  number  of  particles  initially  in  an  auroral  flux  tube  with 

2 

an  area  of  1  cm  in  the  auroral  ionosphere  will  be  defined  here  as 

N  =  Vn  where  n  is  the  particle  density  referred  to  above  and  V  is  the 

"effective  volume"  of  the  flux  tube.  For  a  magnetic  flux  tube  connected 

^  12  3 

with  the  auroral  zone  V  *  10  cm  between  the  ionosphere  and  the 
equatorial  plane.  Now  from  (13)  and  the  definition  of  N, 


dN  =  N  f(n,B  )  dB  dju  (14) 

m  m 

where 

fCp.B,,)  -  (B J2)  (Bn  -  b//2] 

Thus  the  particles  within  the  flux  tube  can  be  described  in  terms  of  a 
distribution  of  magnetic  moments  tod  mirror  fields. 

The  value  of  the  magnetic  field  in  the  auroral  ionosphere  B  can 

ft 

be  taken  as  a  lower  limit  for  trapping.  If  an  electric  field 

exists  along  B  such  that  the  mirror  field  for  a  particle  becomes  greater 

than  B  ,  the  particle  can  be  said  to  be  discharged  into  the  atmosphere. 

ft 

That  is,  from  Eq.  (11)  if  initially  the  mirror  field  B  is  greater  than 

in 

Bm  "  (1  ■  Ba  =  to'/n)  Ba 


for  a  given  \i ^  defining  toy  (10),  the  particle  is  discharged, 
number  of  particles  still  trapped  at  a  given  value  of  E^  is  thus 


toy 


nBa  •‘Vi* 


f(u,B  )  dp  (u’/u)  <JB' 

01  m 


B 

0 


The 

given 


(15) 
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Hie  time  rate  of  change  of  N'  exactly  corresponds  to  the  auroral  flux 
Hius  differentiating  (15),  it  can  be  shown  that 


H  ^  Ba  I  *  f(„.  Hi  Ba)  (16) 

““i 

where  p^,  is  the  rate  of  change  of  the  magnetic  moment  p^  that  defines 
Eg.  Now  the  smallest  magnetic  moment  remaining  in  the  flux  tube  is 
given  by  p..  =  p.  B  /B  -  B  ),  since  the  equatorial  magnetic  field  is  not 
zero.  This  moment  defines  the  lower  limit  for  the  integral  in  Eq.  (l6). 
Hie  integrand  in  (l6)  is  determined  from  Eq.  (14).  Suitable  change  of 
variables  leads  to 


$>  =  N 


g(x)  dx 

2X1/2  (bx -  a)1/2 


(IT) 


where  x  =  (p  -  p^)  B&,  =  P^BQ/(l  -  ®0/Ba)>  B  =  U  -  B0/Ba)>  and 

a  =  p^Bo. 

Now  g(x)  can  be  given  for  a  kinetic-energy  spectrum  of  the  form 
W  where  k  is  an  integer.  In  this  case 


g(x)  =  2(k  -  1)  W w"k 


where  corresponds  to  x^  and  is  given  by  =  a/b  as  defined 

above.  is  the  smallest  particle  kinetic-energy  in  the  distribution. 

Hence, 


"Hf  <k  *  l> 


xk+1/2  (tT-  IF'2 


(18) 


and  on  integration, 
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«  =  N(1  -  B  /B  ) 
x  o'  a' 


^  KOO  ^ 


(19) 


where  K(k)  is  given  by 


K(k) 


(2k  -  2) 

(2k  -  4).  ..(2) 

(2k  -  1) 

(2k  -  3).  ..  (3) 

k 


1) 


(20) 


For  k  =  5,  K(k)  =  1024/315*  while  for  k  =  2,  K(k)  =  2.  The  energy 
spectrum  is  evidently  not  an  overwhelming  factor  in  determining  the 
auroral  particle  flux  produced  by  a  given  particle  density,  electric 
field,  and  time  rate  of  change  of  electric  field. 

In  order  to  determine  the  time  scale  of  an  auroral  flux  and  of  the 
electric  field  E 2  alongJ3,  Eq.  (19)  can  be  written 


N  =  -  Np  (21) 

where  0  =  K(k)  (1  -  BQ/Ba)^  .  Thus,  integrating  (2l),  the  integral  of 
the  flux  <J>  over  time  for  a  given  display  can  be  derived  as 


ft 

1  <b  dt  =  N  -  N  = 

fi  -  (h/m^  )'P] 

'0 

L  0  J 

(22) 


where  Nq  is  the  initial  number  of  particles  in  the  flux  tube,  N  is  the 

number  at  time  t,  is  the  initial  smallest  moment,  is  the  smallest 

o 

moment  at  time  t,  and  (3  is  as  given  above.  Therefore,  the  time  required 

for  an  auroral  flux  of  constant  magnitude  to  deplete  a  region  of 

particles  with  moments  less  than  100  will  be  about  Nq/$.  If 

N  =  KT  cnr  •  10  particle  s/cm J  and  9  =  Kr  particles/cm  sec,  then 
4 

t  *  10  sec,  or  about  3  hours. 
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The  flux  of  atmospheric  particles  up  magnetic  field  lines  is 
neglected  in  the  above  calculation  of  auroral  flux.  This  calculation 
could  be  considered  applicable  to  electron  excesses  in  the  trapped 
plasma.  For  electron  excesses,  the  flux  of  electrons  down  field  lines 
would  not  necessarily  exceed  the  flux  of  atmospheric  protons  up  field 
lines.  The  electron  flux  would,  however,  contribute  the  energy 
associated  with  the  visible  aurora.  Hence  auroral  electron  flux 
estimates  based  on  calculations  from  the  spectra  of  aurora,  or  measured 
directly  in  the  ionesphere,  must  be  taken  as  implying  lower  limits  to 
the  rate  at  which  excess  charge  is  supplied  to  a  given  magnetic  flux 
tube.  For  a  particle  flux  of  10  particle s/crn  sec  in  the  auroral 
zone,  the  lower  limit  to  the  rate  at  which  excess  charge  is  supplied 
near  the  equatorial  plane  (in,  say,  l/lO  the  volume  of  the  flux  tube) 
would  be  about  10  particles/lC;  cnr  sec  or  10_  particles/cm'3  sec. 
This  rate  of  charge  separation  could  be  obtained  with  drift  velocities 
of  the  order  of  100  m/ sec  and  plasma  density  gradients  of  the  order  of 
10  3  ( parti cle/cnr)/km.  If  the  process  of  formation  of  thin  auroral 
sheets  proposed  in  Section  II  is  effective,  the  auroral  flux  would  be 
restricted  to  localized  regions  connected  along  to  the  boundaries  of 
irregularities  in  the  trapped  plasma. 

A  further  result  of  the  electric-field  model  for  discharge  of 
trapped  particles  is  that  some  variation  of  the  energies  of  auroral 
particles  as  a  function  of  time  would  be  anticipated.  Solving  Eq.  (22) 


(23) 
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Now  the  moment  is  given  very  nearly  by  w/b&,  where  W  is  the  discharged 

energy  of  a  particle  and  B  is  the  atmospheric  magnetic  field.  Similarly 

& 

\i ^  =  Wo/Ba  where  Wq  is  the  initial  energy  of  a  discharged  particle  at 

o 

time  t  =  0.  Note  that  Wq  is  not  necessarily  zero.  Thus  at  time  t,  the 
energy  of  a  particle  mirroring  in  the  atmosphere  is  given  by 


W  =  W 

o 


r 

1  -  (l/N  )|  $>  dt 


-i/e 


(24) 


Hie  quantity  in  brackets  goes  to  zero  as  t  becomes  large.  Thus,  as 
particles  of  lower  energy  are  discharged  and  as  the  component  of  the 
electric  field  parallel  to  ^  increases  with  time,  the  energies  of 
particles  producing  a  given  auroral  display  increase  with  time. 
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IV.  POLAR -ELECTROJET  CURRENT  SYSTEMS 


As  mentioned  in  Section  I,  charge  separation  in  the  geomagnet ically 
trapped  plasma  has  been  invoked  to  account  for  the  occurrence  of  polar - 
electrojet  current  systems  (Vestine,  19*50;  Chamberlain,  Kern,  and 
Vestine,  i960;  Kern  and  Vestine,  196I;  Kern,  1961).  Die  process  of 
charge  separation  described  in  Section  II  transports  charge  across 
plasma  density  irregularities  in  the  trapped  radiation.  The  electric 
fields  and  potentials  associated  with  separated  charge  are  short- 
circuited  by  conduction  along  magnetic  field  lines  into  the  atmosphere. 
Urns,  atmospheric  current  systems  would  be  associated  with  such  charge 
transport.  Allowing  magnetic  field  lines  to  be  equipotential  lines  does 
not  change  this  situation.  A  finite  resistance  to  particle  flow  in  the 
atmosphere  ensures  differences  of  potential  between  regions  of  excess 
positive  and  negative  charge  in  the  trapped  radiation.  Die  system  can 
be  visualized  as  involving  a  specified  current  applied  to  a  medium  of 
anisotropic  conductivity;  i.e.,the  atmosphere. 

If  polar-electro Jet  current  systems  are  ascribed  to  Hall  conduction 
caused  by  meridional  electric  fields  in  the  E-region  (Baker  and  Martyn, 
1953;  Chamberlain,  Kern  and  Vestine,  i960;  Akasofu,  private  communication), 
eastward-directed  electrojets  are  associated  with  poleward-directed 
electric  fields  and  westward-directed  electrojets  with  equatorward- 
directed  electric  fields.  Die  drift  directions  produced  by  eastward- 
directed  magnetic-field  gradients  are  such  as  to  produce  poleward- 
directed  electric  fields,  while  westward-directed  magnetic  gradients 
produce  equatorvard-directed  electric  fields  in  the  trapped  radiation 
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incident  in  auroral  regions.  Eastward-directed  electrojets  may  therefore 
be  associated  with  eastward-directed  magnetic-field  gradients,  and 
westward-directed  electro jets  with  westward- directed  gradients  (as  in 
Fig.  3). 

The  extent  in  longitude  of  such  meridional  electric  fields  would 
correspond  to  the  extent  of  the  eastward-  or  westward-directed  magnetic- 
field  gradients  where  irregularities  in  plasma  density  exist  to  allow 
the  development  of  polarization.  Electrojets  would  develop  where 
sufficient  trapped  radiation,  polarized  and  discharged  by  charge 
separation,  penetrates  the  ionosphere.  A  hypothetical  distribution  of 
eastward  and  westward  components  of  the  geomagnetic-field  gradient  in  a 
shell  of  trapped  particles  is  shown  in  Fig.  4.  Adiabatic  drift  of 
electrons  before  midnight  in  this  model  is  northward  to  higher-altitude 
mirror  points  for  a  given  mirror-point  field;  after  midnight,  the 
electron  mirror  points  drift  south  to  lower  altitudes.  TJie  increase  of 
ionospheric  conductivity  encountered  by  the  electrons  on  penetration  to 
lower  altitude  would  seem  to  imply  that  polar  electrojets  due  to  Hall 
conduction  will  have  greater  intensity  after  midnight  than  for  a  positive 
bay  current  system  developing  before  midnight. 

Bie  process  described  above  requires  that  the  direct  conduction 
current  associated  with  the  Hall-current  electrojets  be  supplied  by 
continuous  charge  separation  in  the  trapped  radiation.  This  current  is 

equal  therefore  to  the  charge  transport,  nev^r/  per  radian  of  longitude, 

_2o 

where  n  is  the  plasma  number  density,  e  =  1.6  x  10  emu  is  the 

electronic  charge,  v^  is  the  separation  velocity  for  particles  of 

q 

opposite  sign,  r  *  5  x  10  cm  is  the  average  radius  of  the  region  of 
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drift  separation  conjugate  to  auroral  latitudes,  and  £  =  2  x  10  ^  cm  is 
the  approximate  extent  of  the  drift  region  along  lines  of  force.  The 
total  direct  conduction  current  per  radian  of  longitude  in  the  auroral 
region  associated  with  a  typical  electrojet  can  be  estimated.  If 

Baker  and  Martyn1 s  (1953)  value  for  the  height-integrated  direct 

-8  4 

conductivity  of  10  emu  and  an  estimated  electric  field  of  10  emu 

g 

are  used,  and  if  the  radius  of  the  auroral  zone  is  taken  as  2.5  x  10°  cm, 

4 

this  current  turns  out  to  be  2.5  x  10  emu  per  radian  of  longitude  in 
the  auroral  zone.  Equating  this  to  the  drift  current  given  above,  and 
taking  n  =  10  particles/ cm  ,  permits  calculation  of  the  drift  separation 

4  / 

velocity  required  to  maintain  such  a  current  as  v  =  1.5  x  10  cm/sec. 

We  can  apply  the  equation  given  earlier  for  the  drift  velocity  to 

obtain  an  estimate  of  the  required  geomagnetic  field  distortion.  Since 

the  magnitude  of  the  east-west  component  of  the  magnetic  field  gradient 

can  be  written^  B  =  (l/r)  dB/d$,  from  (l)  dB/d®  =  B2ervdW(l  +  cos2  a). 

With  B  =  10  ^  gauss  at  7  to  8  earth  radii,  e  =  1.6  x  10-20  emu, 

9  4  ,  r 

r  =  5  x  10  cm,  v^  =  1.5  x  10  cm/sec,  W  =  6  Kev,  and  taking  the  average 
value  for  a  of  n/4,  we  obtain  dB/d®~10-1+  gauss/radian.  This  amount  of 
distortion  would  be  considered  applicable  to  times  of  magnetic  bays. 

Note  here  that  the  energy  density  of  the  particles  required  in  the 
trapped  plasma  is  of  the  order  of  60  Kev/cm^.  Ifcis  is  about  1/3  the 
energy  density  of  the  geomagnetic  field  at  7  earth  radii.  If  the  energy 
density  of  the  particles  is  higher,  smaller  distortions  could  lead  to 
the  necessary  charge  transport.  When  the  energy  densities  are  less,  on 
the  other  hand,  the  distortions  of  the  field  must  be  greater  in  order  to 
transport  sufficient  charge  to  drive  the  electrojets. 
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V.  TEE  POLAR  AURORA 

Hie  adiabatic  charge-separation  mechanism  described  in  Section  II, 
together  with  the  discharge  mechanism  discussed  in  Section  III,  provides 
a  means  for  producing  the  polar  auroras  from  a  plasma  trapped  in  the 
geomagnetic  field.  The  results  indicate  that  auroral -particle  fluxes 
will  depend  on  (l)  the  energy  density  of  particles  in  the  plasma,  (2)  the 
distortion  of  the  geomagnetic  field,  and  (3)  plasma  density  gradients 
associated  with  irregularities  in  the  trapped  plasma.  Hie  energies  of 
incident  auroral  electrons  are  found  to  increase  as  a  function  of  time 
because  of  the  dependence  of  auroral -particle  flux  on  the  distribution 
of  particle  energies  and  magnetic  moments  within  the  plasma  connected 
along  B  with  the  auroral  region. 

In  Section  II,  it  was  pointed  out  that  the  formation  of  thin 
auroral  sheets  might  be  a  natural  consequence  of  the  proposed  charge 
separation.  Hiis  geometry  has  been  considered  by  Kern  and  Vestine  (1961), 
who  point  out  some  of  the  consequences  of  charge  separation.  Among  these 
consequences  is  the  fact  that  sheet  beams  of  charged  particles  are 
subject  to  electrodynamic  instabilities  as  demonstrated  both  experimentally 
and  theoretically  (Kyhl  and  Webster,  1956;  Pierce,  1956;  Webster,  1955* 
1957)*  Such  considerations  can  be  applied  to  auroral  morphology  to 
account  for  the  observed  sequences  of  auroral  forms  in  terms  of  changes 
in  the  velocity  and  flux  density  of  incident  auroral  particles. 

Hie  principles  discussed  above  can  be  used  to  interpret  observations 
of  auroras  and  electrojet  current  systems.  For  example,  one  form  of 
distortion  of  the  geomagnetic  field  might  be  such  as  to  produce 
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components  of  the  geomagnetic  field  gradient  directed  toward  the  midnight 
meridian  in  the  equatorial  plane  on  the  night  side  of  the  earth.  Such 
gradients  can  be  thought  of  as  lying  in  surfaces  of  constant  plasma 
density  in  the  magnetosphere.  If  this  is  the  case,  charge  separation 
of  the  kind  discussed  in  Section  II  will  lead  to  polarization  of  the 
trapped  plasma.  Hie  pattern  described  here  is  shown  in  Fig.  4.  If  the 
region  of  trapped  plasma  is  limited  in  extent,  prior  to  midnight  an 
electron  excess  will  accumulate  on  the  outer  boundary.  After  midnight, 
a  proton  excess  will  appear  on  the  outer  boundary. 

Hie  systems  of  currents  and  the  meridional  electric  fields 
associated  with  charge  separation  of  the  kind  shown  in  Fig.  4  are  such 
as  to  account  for  many  features  of  the  magnetic  observations  of  Heppner 
(1954)  for  College,  Alaska.  Observed  drift  of  auroral  structures  are 
accounted  for  in  terms  of  motor  drift,  westward  drift  being 
produced  before  midnight  and  eastward  after  midnight  (Nichols,  1957} 

Kim  and  Currie,  1958).  Hie  drift  of  incident  protons  southward  before 
and  northward  after  midnight,  as  observed  by  Reid  and  Rees  (1961),  is 
also  accounted  for  by  this  configuration. 

Auroral  morphology  and  drift  directions  of  Davis  (1962)  are  shown 
in  Fig.  5.  Hiese  distributions  of  auroral  forms  imply  an  east— west 
charge  separation  derived  from  Stflrmer-orblt  motion  in  the  geomagnetic 
field  combined  with  radial  drift  separation  due  to  eastward  or  westward 
components  of  the  geomagnetic- field  gradient  shown  in  Fig.  4.  Hie  fixed 
spatial  relation  of  the  auroral  pattern  with  respect  to  the  sun  observed 
by  Davis  (1962)  is  of  course  basic  to  the  idea  that  charge  separation  is 
caused  by  interaction  of  the  magnetosphere  and  an  ionized  solar  stream. 
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Heupner  (1954)  shows  that  auroral  forms  coinciding  with  peak 
westward-directed  electro jets  are  often  active  rayed  forms.  Such  forms 
would  be  consistent  with  the  model  developed  above,  in  that  they  would 
coincide  with  peak  electrostatic-field  development  in  the  trapped 
radiation.  Hie  activity  would  be  ascribed  to  electrodynamic  instability 
arising  from  the  characteristic  behavior  of  beams  of  energetic  charged 
particles  (Webster,  1957}  Kern  and  Vestine,  1961).  The  peak  electrojet 
current  would  be  similarly  associated  with  the  highest  electric  field  in 
the  trapped  radiation. 

Irregularities  in  plasma  density  are  required  to  generate  thin 
auroral  structures.  The  assumption  that  such  irregularities  exist  in 
the  neutral  plasma  connected  with  auroral  regions  is  supported  by 
observed  fluctuations  in  the  counting  rates  of  satellites  and  probes 
passing  through  this  region,  which  have  been  interpreted  to  indicate  just 
such  structures  (Rosen,  Farley,  and  Sonett.  i960) .  Direct  evidence  for 
the  "dumping"  of  trapped  radiation  has  been  obtained  through  correlation 
of  Explorer  VII  (1959  iota)  radiation  measurements  with  observations 
of  visible  and  subvisible  auroral  emissions  (O’Brien,  et  al.,  i960) . 
Additional  evidence  for  the  dumping  of  trapped  radiation  has  been  obtained 
from  an  analysis  of  satellite  observations  by  Rothwell  and  Mcllwain 
(i960). 

Electric  fields  parallel  to ^ have  been  invoked  by  McDiarmid,  Rose , 
and  Budzlnskl  (1961)  to  interpret  their  rocket  measurements  of  pitch 
angle  and  energy  distributions  in  the  primary  electron  flux  in  an  aurora. 
Hieir  results  also  indicate  space  and  time  separation  of  electron  and 
proton  fluxea.  This  would  support  the  charge-separation  hypothesis. 
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Observations  of  sunlit  auroras  at  high  altitudes  give  additional  evidence 
for  the  existence  of  electric  fields  parallel  to  B.  Hie  presence  in 
such  auroras  of  ions  with  number  densities  of  the  order  of  10  /cm"5 
is  inferred  from  the  fact  that  they  are  bright  enough  to  be  seen.  Bates 
(i960)  suggests  that  such  ions  must  propagate  upward,  since  they  are 
not  normally  present  in  these  regions  in  such  large  numbers.  Upward 
motion  of  positive  ions  is  to  be  expected  in  auroral  structures  produced 
by  the  electric  field  mechanism  described  earlier. 

Horseshoe  auroral  arcs  can  be  explained  as  a  consequence  of  charge 
separation  occurring  in  irregularities  in  the  trapped  neutral  plasma 
which  are  limited  in  east  —  west  extent.  Numerous  so-called  "horseshoe 
aurora"  have  been  reported.  Stflrmer  (1955)  shows  two  illustrations  of 
this  auroral  fbrtn.  Analyses  of  auroral  observations  have  indicated 
that  this  auroral  form  is  not  rare,  but  is  somewhat  difficult  to  photo¬ 
graph.  Observations  available  from  probe  and  satellite  counting  rates 
indicate  that  irregularities  in  the  plasma  number  density  may  indeed  exist 
(Rosen,  Farley,  and  Sonett,  i960).  Given  such  irregularities,  it  can  be  seen 
that  adiabatic  east— west  particle  drift  will  lead  to  charge  separation. 

Uiis  situation  is  represented  in  an  equatorial  section  in  Fig.  6>  the 
charge  separation  depicted  leads  to  a  charge-density  distribution  that 
is  horseshoe-shaped.  Increasing  electric  fields  will  lead  to  precipitation 
in  a  horseshoe-shaped  arc  in  the  auroral  zone.  If  the  mechanism  for 
auroral  6heet  formation  based  on  polarization  drift  is  effective,  such 
arcs  can  be  very  thin  compared  to  the  plasma  irregularities  (see  Section 
II). 


VI.  RECURRENCE  OF  POLAR  AURORA  AND  ELECTOOJETS 


It  has  been  shown  that  the  polar  aurora  and  electro jets  may  arise 
as  a  product  of  adiabatic  particle  motion  in  a  geomagnetic  field  distorted 
by  a  solar  stream.  Recurrence  of  electrojets  and  aurora  at  about  the 
same  longitude  on  successive  nights  can  be  expected  if  there  are  low- 
energy  particles  in  the  trapped  radiation  contiguous  to  the  auroral 
regions.  Low-energy  particles  drift  very  slowly  in  longitude  with 
respect  to  the  earth's  surface,  and  would  effectively  rotate  with  the 
geomagnetic  field.  Such  a  low-energy  plasma  would  exhibit  the  pattern 
of  adiabatic  drift  separation  of  charge  imposed  by  solar-stream  distortion 
of  the  geomagnetic  field.  Local  sequences  of  polar  aurora  and  electrojet 
phenomena  can  therefore  be  expected  to  reflect  a  24-hour  period  corres¬ 
ponding  to  the  rotation  period  of  the  geomagnetic  field. 

Such  distortion  has  not  yet  been  completely  treated  analytically. 

It  would  seem  that  the  atmospheric  phenomena  of  polar  aurora  and  electro- 
jets  may  provide  some  information  regarding  the  nature  of  the  solar- 
stream  geomagnetic -field  interaction.  For  example,  the  hypothetical 
east— west  geomagnetic-field  gradients  proposed  (see  Fig.  4)  to  account 
for  the  auroral  morphology  patterns  of  Davis  (1962)  indicate  that  solar 
streams  may  increase  the  geomagnetic  field  near  midni^vt.  Such  patterns 
may  also  be  accounted  for  by  the  motion  of  field-line  shells  with  respect 
to  surfaces  of  constant  field. 

If  curvature  of  field  lines  is  considered,  the  pattern  proposed  here 
would  correspond  to  the  warping  of  field  lines  out  of  meridional  planes 
toward  the  midnight  meridian.  Such  warping  might  result  if  solar  streams 
tend  to  extend  lines  of  force  on  the  night  side  of  the  magnetosphere. 
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CHARGE  (ELECTRONS) 

VB:  MAGNETIC  FIELD  GRADIENT,  TRANSVERSE  TO  FIELD 
GIVING  RISE  TO  TRANSVERSE  FORCE  PER  PARTICLE 
f-ftVB 

e:  ELECTRIC  FIELD  of  SEPARATED  positive  and 

NEGATIVE  CHARGE 

VE:  NEUTRAL  £  x  fi  MOTOR  DRIFT  VELOCITY  DUE  TO 
ELECTRIC  FIELD  E 

Fig.  I — Charge  separation  in  a  bounded  plasma  due 
to  transverse  geomagnetic  field  gradient 
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INITIALLY  NEUTRAL  PLASMA 
BOUNDED  BY  GEOMAGNETIC 
FIELD  LINES 


Fig.  2 — Charge  separation  in  geomagneticolly  trapped  radiation 
due  to  eastward  geomagnetic  field  gradient 
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Fig.  3  —  Polarization  of  radiation  incident  in  the  auroral  zone 
and  the  Hall  conduction  polar  electrojet  currents 
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Fig.  4  —  Hypothetical  eastward  and 
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-  =  THIN  AURORAL  SHEETS 

=  DIRECTION  OF  MOVEMENT 


Fig.  5 — Auroral  morphology  of  Davis  (i960) 
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Fig.  6  —  Charge  separation  leading  to  horseshoe  auroral  arcs 
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